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ABSTRACT 

We present a statistical study of the environments of massive galaxies in four redshift bins between 
z = 0.04 and z = 1.6, using data from the Sloan Digital Sky Survey (SDSS) and the NEWFIRM 
Medium Band Survey (NMBS). We measure the projected radial distribution of galaxies in cylinders 
around a constant number density selected sample of massive galaxies and utilize a statistical 
subtraction of contaminating sources. Our analysis shows that massive primary galaxies typically 
live in group halos and are surrounded by 2 to 3 satellites with masses more than one-tenth of the 
primary galaxy mass. The cumulative stellar mass in these satellites roughly equals the mass of 
the primary galaxy itself. We further find that the radial number density profile of galaxies around 
massive primaries has not evolved significantly in either slope or overall normalization in the past 
9.5 Gyr. A simplistic interpretation of this result can be taken as evidence for a lack of mergers in 
the studied groups and as support for a static evolution model of halos containing massive primaries. 
Alternatively, there exists a tight balance between mergers and accretion of new satellites such that 
the overall distribution of galaxies in and around the halo is preserved. The latter interpretation is 
supported by a comparison to a semi-analytic model, which shows a similar constant average satellite 
distribution over the same redshift range. 

Keywords: Galaxies: elliptical and lenticular, cD - Galaxies: groups: general 



1. INTRODUCTION 

The properties of galaxies in the nearby universe are 
strongly affected by the environment in which they re- 
side. Morphology, mass, star formation rates and stellar 
colors of individual galaxies have all been shown to be 
correlated with the local galaxy density (e.g., Oemler 
1974; Dressier 1980; Kauffmann et al. 2004; Hogg et al. 
2004; Blanton et al. 2005; Thomas et al. 2005; Clemens 
et al. 2006; van den Bosch et al. 2008; Skibba & Sheth 
2009; Peng et al. 2010; Zehavi et al. 2011; Quadri et al. 
2012). Therefore, significant effort has been devoted to 
studies of galaxy environments at low and high redshift, 
utilizing a number of approaches (e.g., Zehavi et al. 2002; 
Madore et al. 2004; Mandelbaum et al. 2006; Holden et 
al. 2007 Gavazzi et al. 2007; Bolton et al. 2008; Wake 
et al. 2008, 2011; Patel et al. 2009, 2011; Cacciato et al. 
2012; Knobel et al. 2012; Diener et al. 2012). Neverthe- 
less, the density and distribution of galaxies in groups, 
where nearly all massive galaxies are expected to reside, 
are still poorly known at z > 1. At this redshift, it is dif- 
ficult to uniquely match galaxies to halos except in the 
most overdense regions. 

Recently, a number of studies analyzed the small-scale 
clustering of satellite galaxies around massive primaries 
(the most massive galaxies in their halos) using statistical 
tools. By utilizing large data sets and statistical back- 
ground subtraction, authors were able to extract satel- 
lite galaxy distributions from observational and numeri- 
cal catalogs (e.g., Masjedi et al. 2006, 2008; Watson et al. 
2012; Tal et al. 2012a; Quilis & Trujillo 2012; Budzynski 



et al. 2012; Neirenberg et al. 2012; Jiang et al. 2012; 
Wang & White 2012; Marmol-Queralto et al. 2013; Guo 
et al. 2013). 

Here we follow a similar approach and derive galaxy 
number-density profiles out to z = 1.6 to study the prop- 
erties and evolution of galaxy distributions around mas- 
sive primaries. An important advantage of this approach 
is that it does not rely on accurate redshift measurements 
of the studied galaxies. 

Throughout the paper we adopt the following cosmo- 
logical parameters: = 0.3, = 0.7 and Hq = 70 
km s^i Mpc"^ 

2. SAMPLE SELECTION 

We utilized data from two public surveys for this study; 
from SDSS (York et al. 2000) and from NMBS (van 
Dokkum et al. 2009; Whitaker et al. 2011). Galaxies 
in three redshift bins in the range 0.4 < z < 1.6 were 
selected from NMBS, a deep imaging program with the 
NOAO Extremely Wide-Field Infrared Imager. The to- 
tal imaging area of NMBS spans roughly 0.39 deg^ over 
two fields and the catalogs include tens of thousands of 
galaxies for which excellent photometric redshifts were 
calculated {az/{l+z) ~ 0.02; Brammeret al. 2009). Low 
redshift galaxies were selected at 0.04 < z < 0.07 from 
the MPA-JHU emission line analysis catalogs^ which are 
based on the seventh data release of SDSS (Abazajian et 
al. 2009). 

^ http:/ /www. mpa-garching.mpg.de/SDSS/DR7/ 
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Figure 1. Sample matching of the four redshift bins. The teal, 
yellow and magenta lines represent the cumulative mass functions 
of Marchesini et al. (2009). We calculated galaxy masses at a 
constant cumulative number density of 4xl0~^ Mpc~^ (gray line) 
and derived a mass-redshift relation (inset figure). We then se- 
lected galaxies in log M bins according to the derived relation. 

2.1. Cumulative number density matching 

Studies of galaxy evolution often match their samples 
such that galaxies at all redshifts have similar observed 
properties (e.g., constant stellar mass, luminosity). How- 
ever, the masses and luminosities of galaxies are expected 
to evolve and change with time. Alternatively, samples 
can be matched based on their number density, regardless 
of galaxy properties. Since cumulative number density is 
not expected to change dramatically even with a non- 
negligible number of mergers (e.g., Gao et al. 2004; van 
Dokkum et al. 2010; Papovich et al. 2011; Patel et al. 
2012; Leja et al. submitted), this approach essentially 
follows the evolution of observed properties of a given 
galaxy population. We utilized the latter technique to 
match massive galaxy populations in the two data sets 
in four redshift bins. 

We started by adopting the model fits to the mass func- 
tions from Marchesini et al. (2009) to calculate the cu- 
mulative number density of three redshift samples as a 
function of stellar mass (Figure 1). We then fit a line to 
the mass-redshift relation of galaxies at a constant num- 
ber density of 4xl0~^ Mpc~"^ and found the following 
mass evolution with redshift: 

log(M^/M0) = -0.16z -I- 11.43 (1) 

The best-fit relation implies doubling of stellar mass in 
the last 9.5 Gyr for these galaxies and it is consistent 
with other mass evolution studies (e.g., van Dokkum et 
al. 2010). We selected primary galaxy candidates in 
logM^ bins of width 0.3 dex with an evolving median 
mass according to Equation 1. Galaxies were considered 
to be "primary" if no other, more massive, galaxies were 
found within a projected radius of 500 kpc. Otherwise, 
they were counted as "satellites" of their more massive 



Table 1 

Galaxy properties in the four selected redshift bins. 



Data sot 


z 


V [Mpc^] 


log(M„ed/M0) 


Np 


SDSS 


0.04<z<0.07 


1.6x10^ 


11.42 


360 




0.4<2<0.9 


9.8x105 


11.33 


44 


NMBS 


0.9<2<1.3 


1.4x10** 


11.25 


89 




1.3<z<1.6 


1.2xl0'5 


11.20 


97 



neighbor. The redshift limits for this study were de- 
termined such that the samples from both surveys were 
complete down to one-tenth of the stellar mass of all se- 
lected primaries. Redshift bin sizes were chosen to have 
roughly equal volumes in the three NMBS bins and to 
maximize the number of galaxies from SDSS, whilst re- 
maining complete in stellar mass. 
A summary of the selected samples is given in Table 1. 

3. RADIAL NUMBER DENSITY PROFILES 

Various techniques are commonly used to quantify 
galaxy environments at low and high redshift. Here we 
perform an analysis of the distribution of objects around 
the selected galaxies and subtract the contribution of 
background and foreground sources in a statistical man- 
ner. One of the advantages of this method is that it does 
not rely on a priori assumptions regarding the total mass 
or size of the host dark matter halo. 

3.1. Profile derivation 

We started by selecting all galaxies in the NMBS cat- 
alog that were within Az < 0.2 of the measured red- 
shift of the massive primaries and within Az < 0.005 in 
the SDSS MPA-JHU catalog. From this selection we ex- 
cluded all galaxies whose mass was less than one-tenth 
of their corresponding primary mass. We then derived 
the radial number density galaxy profiles by binning the 
selected galaxies in log-radial bins (blue lines in Figure 

2)- . . 

The resulting radial number density profiles include 

galaxies that are physically associated with each primary 
as well as sources in the background and foreground. In 
order to account for this contamination we repeated the 
derivation in randomly selected positions within the en- 
tire survey area and calculated the average radial profile 
of contaminating sources (red lines in Figure 2). Finally, 
we subtracted the average random profile from the aver- 
age primary profile at each redshift to derive the cylindri- 
cal distribution of physically associated galaxies around 
massive primaries. 

The radial distribution of galaxies around massive 
primaries out to z = 1.6 is shown in Figure 3, where 
the density of physically associated galaxies is plotted 
as a function of projected physical distance from the 
primary. The blue line represents the measurement at 
0.04 < z < 0.07 from SDSS while the yellow, green and 
red lines show the three higher redshift bins from NMBS. 
We note that the lowest redshift profile (blue line) is not 
plotted at r < 55" (roughly 74 kpc at z = 0.07), where 
the spectroscopic catalog of SDSS is incomplete due to 
fiber collisions. 
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Figure 2. Demonstration of the radial number density profile extraction technique. We divided catalog sources into log(r) bins around 
the selected primaries (blue lines) and around randomly selected positions in the field (red lines). The bottom-right panel shows an example 
of the stacked radial profiles of all galaxies at 1.3 < z < 1.6. The projected radial profiles of galaxies around the selected primaries were 
revealed by subtracting the average random profile from the galaxy centered profile. 



3.2. Error estimates 

The range of measured number density values in the 
random profiles reflects the statistical variation of the 
overall galaxy distribution throughout the survey fields. 
In addition to averaging all random profiles in each red- 
shift bin, we also calculated the standard deviation of the 
distribution of individual profiles as a function of radius. 
This measurement depicts the statistical uncertainty in 
our calculations. 

In addition to statistical errors, a number of system- 
atic uncertainties may potentially affect our results in a 
non-trivial way. For example, uncertainties in mass and 
redshift measurements may move galaxies in and out of 
a given redshift sample and therefore change the profile 
normalization and slope. In the innermost radii the pro- 
files may be underestimated due to source blending close 
to the bright primaries. While redshift errors in NMBS 
are well quantified and are expected to have a minor 
effect on the number density profiles, uncertainties in 
stellar mass measurements are complex and depend on 
a large number of fitting parameters (e.g., Marchesini et 
al. 2009; Muzzin et al. 2009; Conroy et al. 2009). 

We quantified the effect of redshift errors by scattering 
all measured redshifts in the NMBS catalog by a normal 
distribution with width ^^/(I -I- z) = 0.02. This value 
represents the scatter in the relation between photomet- 
ric and spectroscopic redshift measurements in NMBS as 
was found by Brammer et al. (2009). We extracted num- 



ber density profiles using the scattered catalog, following 
the procedure described in Section 3.1. The overall shape 
and normalization of the resulting profiles are consistent 
with the non-scattered profiles and the resulting errors 
are within ±20% of the above derived statistical errors at 
all radii. The error bars in Figure 3 represent statistical 
uncertainties as well as redshift errors in this study. 

Measured values and error estimates that are presented 
in Figure 3 are given in Table 2. 

3.3. Profile dependence on primary mass 

In order to test whether our results are sensitive to 
the mass of the primaries we divided the 0.9 < z < 1.3 
sample into four bins according to the total stellar mass 
of the primary galaxies. We selected primaries with 
masses in the range 10.85 < log (Mp/M©) < 11.25 and 
calculated the radial number density profile of each mass 
bin following the method described in Section 3. We 
estimated the statistical uncertainties for each mass bin 
using the variation in individual profile measurements. 
Figure 4 shows that the radial number density profiles 
of all four mass bins are consistent with one another and 
no strong mass dependence is apparent in this sample. 
However, we note that the relatively small mass range 
of selected primaries may be insufficient to detect even 
a moderate mass dependence. 



3 



The Astrophysical Journal, 000:000 



Tal et al. 



ID- 



IO- 



CY 



10- 



10- 



10- 




0.04 < z < 0.07 
0.4 < z < 0, 
0.9 < z < 1.3 
1.3 < z < 1.6 
GuoetaL (2011) 
Range of virial radii 



10 



100 



1000 



Radius [kpc] 



Figure 3. Projected radial profile of galaxies around the selected 
primaries in four redshift bins (colored lines). The profiles show a 
remarkable lack of evolution in the radial galaxy distribution since 
z = 1.6. The shaded gray profile represents a similar derivation 
using the semi-analytic model of Guo et al. (2011), over the same 
redshift range. The vertical light gray area shows the range of 
virial radius values of the dark matter halos containing the modeled 
galaxies. The observed profiles are consistent with the modeled 
ones over most of the radial range. 



3.4. Environment of massive primaries 

In order to calculate the mass density in the vicinity 
of our sample primaries we converted the number den- 
sity profiles from Figure 3 into a cumulative number of 
satellites as a function of radius. Figure 5 shows the in- 
tegrated curves of the NMBS profiles as a function of 
integration radius. The top x-axis of the Figure shows 
the orbital timescale of a test particle in a spherical halo 
with total mass 2 x 10^"^ Mq which is distributed as a 
Navarro, Frenk & White (1997; NFW) profile. The y-axis 
on the right-hand side of Figure 5 represents an estimate 
of the total enclosed mass as a function of radius, assum- 
ing that the mass of each plotted galaxy equals 40% of 
its primary mass. This value represents the median mass 
of all selected galaxies as calculated from the global mass 
function sX z = 1.6. 

Average properties of massive galaxy environments out 
to z = 1.6 can be directly measured from these integrated 
profiles. For example, these profiles can be directly com- 
pared with estimates of galaxy pair fractions (e.g., Le 
Fevre et al. 2000; Patton et al. 2002; Lin et al. 2004; 
Belletal. 2006; Kartaltepe et al. 2007; Blucketal. 2009; 
Bundy et al. 2009; Lopez-Sanjuan et al. 2012). Recently, 
WiUiams et al. (2011) and Newman et al. (2012) found 
that roughly 20% of all massive galaxies have a luminous 
companion with a mass ratio of up to 10:1 and within a 
projected distance of 30 kpc. The same result can be di- 
rectly extracted from Figure 5 by reading off an average 
value of 0.2 galaxies within 30 kpc of the primaries. 

In addition to quantifying pair fractions in discrete sep- 
aration values, the derived cumulative profiles describe 



10.85 < log(Mp/Mo) < 10.95 
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Figure 4. Projected radial number density profiles in the redshift 
range 0.9 < z < 1.3 in four mass bins. We found no evidence for a 
strong dependence on the primary galaxy mass within the probed 
mass range. 

the enclosed mass in a range of radii, essentially quanti- 
fying the average stellar mass content of studied halos. 
The total number of galaxies within a mass range of 1:10 
and within roughly 400 kpc of the primary is on average 2 
to 3 in all redshift bins. This value suggests that massive 
primary galaxies have typically resided in group halos in 
the last 9.5 Gyr and it is comparable to the number of 
luminous satellites in massive groups from other studies 
(e.g., Tal et al. 2012a; Quilis & Trujillo 2012), Moreover, 
the total stellar mass in satellites within the halo virial 
radius can be inferred from the cumulative profiles and 
it roughly equals the mass of the primary galaxy itself. 

3.5. Profile redshift evolution 

The average radial profiles of the four redshift bins 
show a remarkable lack of evolution in shape and overall 
normalization. Within the measurement scatter, the pro- 
files are consistent with one another on nearly all scales 
from z = 1.6 to z ~ 0. The lack of evolution in the 
number density profiles can be naively taken as evidence 
for a static model of massive halo evolution. According 
to the stable clustering model (Davis & Peebles 1977), 
cluster halos decouple from the large-scale cosmological 
flow shortly after formation. By that point, the halos 
are virialized and they show no significant evidence for 
radial infall. As a consequence, any mergers that involve 
primary galaxies can only be with other galaxies from 
the same halo. 

In order to test the effect of in-halo mergers on the 
evolution of the radial number density profile we sim- 
ulated mergers starting with the highest redshift bin 
(1.3 < z < 1.6) and compared the results to the mea- 
sured profiles at lower redshifts. Mergers were simulated 
by excluding galaxies from the radial galaxy distribu- 
tions in the range 10 < r/kpc < 400, essentially allowing 
satellites to merge with the primary from anywhere in 
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Figure 5. Cumulative number of galaxies as a function of inte- 
gration radius from NMBS (solid yellow, green and red lines) and 
from Gil (gray shaded region). The upper x-axis shows an esti- 
mate of the orbital time scale under an assumed NFW halo with 
total mass 2 X 10^^ Mq. The right-hand y-axis displays the ra- 
dially integrated galaxy mass, assuming that all galaxies have the 
median mass of 40% of their respective primary mass. Primary 
galaxies can double in mass by merging with 2.5 galaxies that are 
at r < 300 kpc. 



the halo. For an estimate of the average merger rate 
we utiHzed the results from van Dokkum et al. (2010), 
who found that massive primary galaxies have roughly 
doubled in stellar mass since z = 2. Assuming that 
most of the mass growth takes place via minor merg- 
ers with mass ratios less than 10:1 (Tal et al. 2012b), 
the implied merger rate is roughly 0.5 Gyr~^. Figure 
6 shows the modeled profiles at the two lowest NMBS 
redshifts (dashed lines), as well as the observed profiles 
(solid lines). 

It is clear from Figure 6 that in-halo mergers would 
have a pronounced effect on the evolution of the radial 
number density profile. If each primary galaxy under- 
went mergers at the estimated rate and no new galaxies 
were added to the halos on similar scales, the observed 
profile at 1.3 < z < 1.6 (solid red line) would evolve into 
the modeled profiles at 0.9 < z < 1.3 (dashed green line) 
and at 0.4 < z < 0.9 (dashed yellow line). Instead, the 
observed profile at these redshifts (solid green and yellow 
lines) are consistent with no evolution at all. 

Taken at face value, this analysis supports a static 
model, suggesting that a significant rate of galaxy merg- 
ers within the halo is unlikely on time scales of up to 
roughly 10 Gyr. However, studies of the stellar mass 
function of massive galaxies show that such galaxies grow 
significantly in mass over the same redshift range, mainly 
through mergers (e.g., Dickinson et al. 2003; Bundy et 
al. 2006; Drory et al. 2005; Pozzetti et al. 2007; Perez- 
Gonzalez et al. 2008; Marchesini et al. 2009). It is 
therefore evident that truly non-evolving models insuffi- 
ciently recover the observed evolution of massive primary 
galaxies and their environments. 
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Figure 6. Estimates of the effect of in-halo mergers on the radial 
number density profile evolution. We modeled mergers by exclud- 
ing galaxies from the derived profiles of each primary galaxy in the 
range 10 < r/kpc < 400. The dashed lines represent the expected 
evolution of the 1.3 < z < 1.6 profile into subsequent redshift bins 
at a merger rate of 0.5 Gyr~^ with no external halo accretion. 

4. COMPARISON TO A SEMI-ANALYTIC MODEL 
4.1. Density profiles and mass growth 

In order to study the evolution of massive galaxy 
groups in the last 9 Gyr we examined numerical predic- 
tions using the semi-analytic model of Guo et al. (2011; 
Gil). Gil applied a semi-analytic model to the merger 
trees of the JVIillennium Simulation (Springel et al. 2005) 
and improved the treatment of gas dynamics and tidal 
disruptions compared to previous studies. 

For consistency, we first test whether our cumulative 
number density matching technique indeed identifies pro- 
genitors and descendants in Gil. We selected modeled 
galaxies at z — 1.63 using our cumulative number den- 
sity selection criteria and followed their stellar mass evo- 
lution at all simulated time points down to z = 0. The 
median growth factor of massive primary galaxies in Gil 
is roughly 2.4, consistent with the growth that was mea- 
sured observationally for similarly massive galaxies by 
van Dokkum et al. (2010). We also followed the simu- 
lated mass evolution of the group dark matter halos and 
compared it to the average ratio between the halo mass 
of galaxies selected at z = and at z = 1.63. We found 
that both simulated and predicted mass growth factors 
are within 4% of one another with an average halo growth 
by a factor of roughly 3.7. In short, by matching mas- 
sive galaxies by their cumulative number density we were 
able to correctly trace progenitor and descendant halos. 

From the galaxy catalogs of Gil we selected all pri- 
mary galaxies at redshifts z =0.06, 0.62, 1.17 and 1.63, 
following the same selection criteria as those used for 
the observed data. We then collapsed the simulation 
box along one of the physical dimensions and extracted 
the projected radial distributions of physically associ- 
ated galaxies around the selected primaries following the 
method described in Section 2. This approach allows 
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Figure 7. Fate of galaxies around massive primaries in Gil. Solid 
lines represent the fraction of modeled galaxies that merge with 
their respective primary by z ~ as a function of initial separation. 
Dotted lines show the derived radial component of the velocity 
vector at the plotted redshift. Three initial redshift measurements, 
2 = 1.6, 1.0 and 0.4 are plotted in red, blue and green, respectively. 
Galaxies outside of 1000 kpc initially (roughly two virial radii) 
merge with the primaries at a low rate of less than 20% and move at 
a high radial velocity. Galaxies inside the halo virial radius merge 
with the primary at a significantly higher rate with an average 
infall velocity that is close to zero. 

for a direct comparison with the observational results, 
as the uncertainties associated with the profile extrac- 
tion procedure are preserved. The gray profile in Fig- 
ure 3 represents this measurement in the full redshift 
range, 0.06 < z < 1.63. The vertical light-gray region in 
the Figure shows the range of virial radius values of the 
dark-matter halos in which all selected primaries reside, 
as determined by Gil. Finally, the gray curve in Figure 
5 shows the integrated number density profiles over the 
same redshift range. It is evident from Figure 3 that the 
average modeled and observed profiles agree well with 
each other on all but the smallest scales. We note that 
the innermost observational data point (at r<10 kpc) is 
likely underestimated since galaxies inside this radius are 
within the full-width half maximum of the average point 
spread function at z > 1. 

4.2. Fate of modeled satellite galaxies 

The good agreement between our result and the theo- 
retical prediction prompted us to further investigate the 
fate of galaxies around massive primaries by following 
their merging histories and orbits in the models. To do 
so, we selected 1000 primary galaxies at z — 1.63 from 
Gil at random, following the selection criteria from Sec- 
tion 2. We then traced the positions of galaxies in 3 
Mpc apertures around the primaries, that were more 
massive than one-tenth of the primary mass. We re- 
peated the selection at all output redshifts in the range 
0.06 < z < 1.63 and followed the trajectories of individ- 
ual galaxies through time. In addition, we determined 
whether a satellite merges with its primary by z = 0.06 
using the merger trees of the Millennium Simulation. 



The solid lines in Figure 7 show the fraction of satel- 
lite galaxies that eventually merge with their primary 
galaxy as a function of physical separation at z = 1.6, 
1.0 and 0.4. Galaxies that are found inside the virial ra- 
dius at any given redshift are significantly more likely to 
merge with their primary ( f merge > 60% dX z = 1.6) than 
galaxies that start off far from the halo {f merge < 20%). 
However, as is shown in Figure 5, only a small fraction 
of all galaxies within 3 Mpc are found inside of the group 
virial radius (~10%). The y-axis on the right-hand side 
of Figure 7 and the dashed lines show the median ra- 
dial velocity of satellite galaxies with respect to their 
primaries, as a function of physical separation. At each 
redshift we estimated the average streaming velocities by 
measuring the mean distance between satellites and pri- 
maries at two consecutive output redshifts and dividing 
the difference between them by the total modeled time 
span. This measurement gives a rough estimates of the 
average radial component of the velocity vector in the 
rest frame of the primary galaxy. Figure 7 shows that 
galaxies outside of roughly two virial radii fall into the 
halo with high velocities. 

The fate of satellite galaxies in the semi-analytic model 
of Gil seems to be bimodal. The vast majority of galax- 
ies fall toward the most massive halo galaxy with a high 
average radial velocity (in agreement with e.g., Tormen 
1997; Benson et al. 2005; Kochfar & Burkert 2006; Wet- 
zel 2011) and do not merge with the primary by z ~ 0. A 
small fraction of galaxies have a significantly lower mean 
radial velocity and they eventually merge with the pri- 
mary. The probability that a galaxy will merge with its 
massive primary is inversely proportional to the physical 
distance between them. At any simulated redshift, most 
galaxies that are found outside of the halo virial radius 
will not merge with the primary while a significant frac- 
tion of the galaxies found inside of the virial radius will 
do so by z ~ 0. 

We note that despite the good agreement between ob- 
served and modeled galaxy profiles, the fate of satel- 
lite galaxies in numerical simulations should be analyzed 
with caution. In the context of this study, the disagree- 
ment between the evolution of modeled and observed 
stellar mass functions (e.g., Fontanot 2009; Kajisawa 
2009; Cirasuolo 2010; Lu 2012; Mutch 2013) suggests 
that galaxy merging may be misrepresented in numeri- 
cal calculations. Therefore, galaxy merging trajectories 
in observed group halos may be significantly different 
from modeled ones. 

5. SUMMARY AND CONCLUSIONS 

The evolution of galaxy properties is determined in 
large part by the environments in which they reside and 
the mass density around them. Here we quantified the 
radial number density profile of galaxies around massive 
primaries in the redshift range 0.04 < z < 1.6. We 
showed that massive galaxies have typically resided in 
group environments in the past 9.5 Gyr and that they 
are on average surrounded by 2 to 3 massive satellite 
galaxies. We also estimated that the cumulative stellar 
mass of halo satellites within the virial radius accounts 
for roughly as much mass as in the massive primary itself, 
suggesting that the potential for mass growth through 
in-halo mergers is significant. 

Mergers within massive galaxy groups, even at a low 
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rate, are expected to dramatically influence the evolu- 
tion of observed galaxy distributions in such halos, un- 
less an influx of galaxies continuously repopulates the 
halo. We compared the number density profiles of mas- 
sive satellite galaxies in four redshift bins and showed 
that the profiles are consistent with no evolution out to 
z — 1.6. Although this result may be naively interpreted 
as evidence for the insignificance of galaxy mergers in 
group halos, the observed mass growth of massive quies- 
cent galaxies over the same redshift implies that mergers 
occur at a non-negligible rate. Moreover, results from 
semi-analytic models suggest that massive primaries con- 
tinuously merge with satellite galaxies in their halo. 

It would also be interesting to compare the observed 
satellite distributions to predictions from hydrodynam- 
ical simulations, such as those of Naab et al. (2007). 
These models successfully reproduce the build-up of the 
outer envelopes of massive galaxies through minor merg- 
ers (e.g., Naab et al. 2009; Oser et al. 2012; Gabor & 
Dave 2012; Hilz et al. 2013), and an important ques- 
tion is whether the reservoir of satellite galaxies is also 
reproduced. 

In conclusion, the observed lack of evolution in the 
number density profiles suggests that there exists a tight 
balance between mergers and accretion in massive galaxy 
halos. As satellite galaxies merge with their massive pri- 
mary, other galaxies get accreted into the halo at a simi- 
lar rate but on more extreme trajectories. The two com- 
peting processes result in a remarkably balanced galaxy 
distribution out to r = 3000 kpc, over the redshift range 
< z < 1.6. 

This is the first time that an analysis of the evolution 
of galaxy number density profiles is performed at this 
redshift and it provides a new observational insight on 
galaxies in group halos. 
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Table 2 

Projected radial profile measurements 
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z < 1.3 
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z < 1.6 


log(R)=' 


log(0)b 




log(0)b 




log(0)b 




log(0)b 




1.2 
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^Central value of the logarithmic radius bin 
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''Logarithm of the estimated error of the average projected galaxy density 
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